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Abstract 

We discuss the impact of near detectors at a neutrino factory both on standard oscillation 
and non-standard interaction measurements. Our systematics treatment includes cross sec- 
tion errors, flux errors, and background uncertainties, and our near detector fluxes include 
the geometry of the neutrino source and the detector. Instead of a specific detector concept, 
we introduce qualitatively different classes of near detectors with different characteristics, 
such as near detectors catching the whole neutrino flux (near detector limit) versus near 
detectors observing a spectrum similar to that of the far detector (far detector limit). We 
include the low energy neutrino factory in the discussion. We illustrate for which measure- 
ments near detectors are required, discuss how many are needed, and what the role of the 
flux monitoring is. For instance, we demonstrate that near detectors are mandatory for 
the leading atmospheric parameter measurements if the neutrino factory has only one base- 
line, whereas systematical errors partially cancel if the neutrino factory complex includes 
the magic baseline. Finally, near detectors with v T detection are shown to be useful for 
non-standard interactions. 
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1 Introduction 



In neutrino physics, three-flavor oscillations have been successful to explain all relevant 
neutrino data, see, e.g., Ref. [1]. In particular, the solar and atmospheric oscillation param- 
eters have been measured with very high precisions, and the reactor mixing angle #i 3 has 
been strongly constrained. Future long-baseline and reactor neutrino experiments will test 
this small angle further, and be sensitive to leptonic CP violation and the neutrino mass 
hierarchy (see Ref. [2] and references therein). The ultimate high precision instrument for 
these purposes might be a neutrino factory [3-6]. Using different baselines and oscillation 
channels, it can basically disentangle all of the remaining oscillation parameters [7-10] in 
spite of the presence of intrinsic correlations and degeneracies [6,11-13]. Furthermore, a 
neutrino factory and other future neutrino oscillation experiments will be sensitive to new 
physics, such as so-called non-standard interactions; see, e.g., Refs. [14-38]. 

The design of a neutrino factory has been put forward and discussed in international stud- 
ies, such as in Refs. [2,39-41]. Especially the most recent study, the International Neutrino 
Factory and Superbeam Scoping Study [2,42,43], has laid the foundations for the currently 
ongoing Design Study for the Neutrino Factory (IDS-NF) [44]. This initiative from about 
2007 to 2012 is aiming to present a design report, schedule, cost estimate, and risk assess- 
ment for a neutrino factory. It defines a baseline setup of a high energy neutrino factory 
with two baselines Li ~ 4 000 km and L 2 — 7 500 km (the "magic" baseline) operated by 
two racetrack-shaped storage rings, where the muon energy is 25 GeV (for optimization 
questions, see Refs. [10,45]). There are no near detector specifications yet, and the system- 
atics treatment is done in an effective way by signal and background normalization errors 
uncorrelated among all channels and detectors. Therefore, there are a number of questions 
currently raised within the IDS-NF: 

• Study of the potential of near detectors to cancel systematical errors. 

• Study of the characteristics of the near detectors, such as technology, number, etc.. 

• Study of the use of the near detectors for searches of new physics. 

In this work, we address several of these questions using a more refined systematics treat- 
ment including cross section errors, flux errors, and background uncertainties. Similar sys- 
tematics studies using the pull method for physics studied have so far been especially used 
in the context of reactor experiments, see, e.g., Refs. [46,47]. A systematics discussion of 
superbeams can be found in Ref. [48] . We also include a low energy version of the neutrino 
factory in the discussion [49-51], which has been recently proposed to reduce the accelerator 
cost in the case of large 9 13 . 

For the near detectors at a neutrino factory, a summary of options can be found in Ref. [42]. 
The near detectors are, for example, supposed to measure the neutrino flux, the neutrino 
beam angle and its divergence, the neutrino energy, and the neutrino cross sections. For the 
technology, liquid argon time projection chambers, conventional scintillators, scintillating 
fiber trackers, gas time projection chambers, silicon detectors, and emulsion detectors are 
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mentioned as possible options. We will not specify a particular detector technology, but in- 
stead postulate that the near detectors can at least measure the (or u^) event rates with 
(at least) the same energy resolution as the far detectors. Since close to the neutrino source 
the geometry of both the source and the detector become important, we introduce qualita- 
tively different classes of near detectors with different characteristics, such as near detectors 
catching the whole neutrino flux (near detector limit) versus near detectors observing a 
spectrum similar to the far detector (far detector limit). We do not include any additional 
properties, such as beam angle measurements, in our treatment. The main purpose of this 
study is not an accurate description of the near detectors with all their properties, but an 
estimate of when near detectors are important from the physics point of view. 

The size of these multi-purpose detectors will typically be small, such as 0(100) kg, and 
the charged current event rates extremely high compared to that in the far detector. To 
a first approximation, the near detectors measure the product of flux and cross sections. 
However, using different types of charged current interactions (such as the purely leptonic 
inverse muon decays with a different threshold of about 10.9 GeV), one can, in principle, 
disentangle the fluxes and cross sections in the near detector. Since these interactions have 
relatively small cross sections, larger near detectors may be favorable at the end. We do 
not enter this level of detail, but only assume that the near detectors measure the inclusive 
charged current cross sections. We choose different sizes for the near detectors, from small 
(200 kg), over "typical" (of the size of the SciBOONE, MINERz/A, NOMAD, or the MINOS 
near detectors) to a large (hypothetical) detector to capture the whole beam. For the flux 
knowledge, we make rather conservative assumptions. For example, the total number of 
muons circulating in a storage ring may be inferred from a Beam Current Transformer, such 
that the flux is known to the level of 10~ 3 . However, if the muon beam divergence is too 
large (and not monitored), additional uncertainties enter the flux. Therefore, we start with 
a much more conservative assumption (2.5%), and illustrate the improved flux knowledge 
from different interaction types in the near detectors and beam monitoring separately. 

Note that the purpose if this work is to provide information for the development of near 
detectors, in particular, for which physics measurements near detectors are required. Of 
course, systematics is closely connected to that. Our study is to be interpreted as a guideline 
for the minimal requirements to near detectors and systematics in order to obtain to the 
physics wanted. 

This study is organized as follows: We discuss the neutrino factory flux in Sec. 2. Then we 
introduce our near detectors in Sec. 3, and compute the corresponding spectra. In Sec. 4 we 
describe our systematics treatment and simulation methods, where details can be found in 
Appendix A. In Sec. 5 we show the physics results for the high energy neutrino factory and 
in Sec. 6 for the low energy neutrino factory, both for standard oscillations. Furthermore, 
we discuss in Sec. 7 the non-standard interactions in the context of the high energy neutrino 
factory, before we summarize in Sec. 8. 
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Figure 1: The unpolarized differential neutrino spectra in the laboratory frame for (left panel) and 
f e (right panel), where — 25 GeV and different off-axis angles are used (see legends, in degrees). 

2 Neutrino factory flux 

The (double) differential decay rates of an unpolarized muon in the laboratory frame are 
given by (see, e.g., Refs. [3-6]) 

Pcoad)^ [377V - 4 T £„„(1 - /3cos0)] , (1) 

P cos 0)E 2 Vb K - 2 1 E^(l - /3cos0)] . (2) 

Here v a (a G {e, /x}) stands for both v a and v a , 7 = E^/m^ =1/ a/1 — /3 2 is the boost factor, 
and is the muon rest mass. The angle 9 is the angle between the travel direction of the 
muon and the observer (in the laboratory frame, at the decay point), which we henceforth 
call off-axis angle. Note that in a racetrack-shaped storage ring, any muon polarization 
averages out with a high precision [42] , which means that we can use the above unpolarized 
spectra. In addition, note that the theoretical knowledge on these fluxes (the calculated 
spectrum can be trusted at the level to 1CT 3 ) is expected to be more precise than needed 
for this study. 

Very importantly, the off-axis angle has to satisfy the flavor-independent constraint 

9<9 cnt with 6 cut (E„)~-J^-l, (3) 

7 V E u 

which is a relationship depending on the neutrino energy. Conversely, the inverse of this 
relationship limits the neutrino energy observable at a certain off-axis angle, i.e., the flux is 



d 2 T Gpirifj, 
dE~dcos9 ~ 24tt 3 

dE„d cos 9 47r 3 
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Figure 2: The unpolarized single differential neutrino spectrum integrated over energy as a function of 
the off-axis angle 8 (in radians) for different values of E^. 

zero for E v > E UjCUt (8) ~ E^/(l + (76 1 ) 2 ), which means that the off-axis angle 9 determines 
the maximal available neutrino energy. It needs to be included in any integration over 
off-axis angle (such as over the detector's surface area) or the energy. We illustrate the 
relationship in Fig. 1, where especially the spectrum faces a sharp cutoff depending on 
the off-axis angle. At E u > E UtCnt (9), the flux is zero. From Eq. (3), for E v = E^, we 
have 9 cut — 0, which means that the high energy part of the beam is going in the forward 
direction. For E u = E^/2, we have 9 cnt ~ 1/7, which is often referred to as the opening 
angle of the beam (we will use a different, energy independent definition below). For small 
neutrino energies, the beam obviously becomes wider than that for high energies, which 
corresponds to the relative enhancement of the off-axis fluxes at low energies in Fig. 1. 

For the following discussion, it turns out to integrate out the energy in order to obtain the 
single differential decay rate (observing Eq. (3)) 



dY 



dcosO 384vr 3 [7(1 - (5 cos < 



(4) 



which we show in Fig. 2. Obviously, the energy-integrated decay rate always peaks at 
9 = 0, and decreases monotoneously with the off-axis angle. We define the (neutrino) 
beam divergence 9 by the equation 



I 

fo" 



dr 



dcos9 



sin 9d9 



7 2 (1 — (3 cos( 



sin 9d9 = A , 



(5) 



where T is the total decay rate (integrated over angle and energy) T 



G 2 F ml 

1927T 3 



This means 



that the fraction A of the total flux is contained in the cone limited by the angle 9. The 
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Beam diverg 


;ence 9 


Beam openin 


g angle 9 






A = 0.90 A 


= 0.99 


A = 0.90 


A = 0.99 




= 25 GeV 


0.0127 


0.0420 


0.000983 


0.000300 




= 4.12 GeV 


0.0769 


0.2538 


0.005966 


0.001821 



Table 1: The beam divergence 9 (as denned in Eq. (5)) and the beam opening angle 9 (as denned in 
Eq. (6)) for two different muon energies and two different fractions A in radians. These numbers are obtained 
from Fig. 2. 

beam divergence can be used to define a near detector which captures the total flux of 
the beam. Note that we imply that the muon beam divergence be much smaller than the 
neutrino beam divergence, since otherwise a sizable amount of events will be lost [42]. 

Similarly, we define the beam opening angle 9 as 

dT , dT 

-e ~ ^ ' "77777)1^=0 , (6) 



dcos9 ,y ~ y dcos9' 

which quantifies the angle over which the flux stays almost constant. It can be used to 
define a near detector which observes a flux similar to that of a far detector, for which the 
flux at 9 = is used. Note that this particular definition of the opening angle does not 
depend on the neutrino energy, which we have integrated out in Eq. (4), whereas Eq. (3) 
depends on the neutrino energy. We show 9 and 9 in Table 1 for two different muon energies. 
Obviously, the larger A becomes, the smaller is 9, and the larger is 9. The actual diameter 
of the beam for these two quantities is finally obtained as 

D ~ 2 x L x 9 (7) 

where L (baseline) is the distance between production point and observer. 

3 Definition and spectra of the near detectors 

Here we define our near detectors and describe the computations of the spectra. 
3.1 Near detector definitions 

As indicated in the introduction, we are interested in the conceptually different detectors 
from the flux point of view, where the sole purpose is the flux x cross section measurement. 
We use the racetrack-shaped storage ring geometry from IDS-NF [44] designed for a high 
energy neutrino factory. This geometry together with possible near detector locations (to 
be defined later) is illustrated in Fig. 3. In general, for a near detector, the far distance and 
point source approximations 

d 2 T d 2 T 

'cos#) ~ — — — 1 6» =0 , L > s (Size of source) (8) 



dE v dcos9 dE u dcos9 

do not hold, and the geometry cannot be neglected. 
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Figure 3: Geometry of the muon storage ring and possible near detector (ND) locations (not to scale). 
The baseline L is the distance between production point and near detector, i.e., d < L < d + s. 

Our near detector definitions will be based on the following assumptions or observations: 

• We operate the near detectors on-axis, because for an off-axis operation, the event 
rates at high energies will be suppressed according to Eqs. (1) and (2). 

• The only purpose of the near detectors is the (from fi~ decays) and the (from 
fi + decays) event rate measurement using the inclusive charged current cross sections. 
Since the muons and anti-muons are assumed to circulate in different directions in the 
storage ring (c/., Fig. 3), two near detectors are required. 

• We do not measure the v e and v e event rates, since we do not need the corresponding 
cross sections in the neutrino factory far detector. 1 Therefore, we avoid to complicate 
the cross section discussion (such as to relate the muon neutrino and electron neutrino 
cross sections). 

• We use the same characteristics as in the far detectors, such as energy resolution and 
binning, for the sake of simplicity. 

• We do not extrapolate the backgrounds from the near to the far detectors, but instead 
use relatively large background uncertainties uncorrelated among all channels. 

• We assume the fiducial volume to be cylindrical. In addition, we assume that the 
detector only sees muon decays from the decay straight (denoted by s in Fig. 3). 2 

• We do not include systematical errors potentially uncorrelated among the detectors, 
such as fiducial volume errors and energy calibration errors. These additional system- 

1 In fact, not even the dominating backgrounds (charge mis-identification and neutral currents) depend 
on the electron neutrino cross sections. 

2 The contributions from the curved sections can be estimated to be at the few percent level and can be 
easily computed once the final geometry is known. 
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Parameter 


ND1 


ND2 


ND3 


ND4 


ND5 


Diameter D 


17 m 


4 m 


4 m 


0.32 m 


6.8 m 


Distance d 


80 m 


1000 m 


80 m 


80 m 


1000 m 


Mass 


450 t 


25 t 


25 t 


0.2 t 


2000 t 



Table 2: Definition of our near detector fiducial (active) volumes in terms of diameter, distance d to the 
end of the decay straight (c/., Fig. 3), and mass. The fiducial volumes are assumed to be cylindrical. If the 
density is about 1 g • cm" 3 (such as for a liquid scintillator), the active detectors will be about 2 m long for 
ND1 to ND4. ND5 is assumed to be OPERA-sized (with a cylindrical shape for the sake of simplicity). For 
the low energy neutrino factory, only ND3 is used. ND4 is a down-scale (tabletop) version of ND2 with the 
same ratio D/L. The ND mass of 25 t is at the upper limit of currently used near detectors. 

atical errors may affect the final results, but they cannot be improved upon by near 
detectors. 

Note that we do not necessarily require a magnetic field or charge identification since our 
distance is too short for the neutrinos to oscillate, but we need sufficient flavor identification 
at the level of (at least) the flux normalization error. Alternatively, the different flavors can 
be distinguished by the different charges, but this may be technically more challenging. 

From the conceptual point of view, we postulate that the geometry determines the limiting 

cases: 

Near detector limit In this case, the neutrino beam divergence given by Eq. (5), applied 
to Eq. (7), is smaller than the detector diameter for the farthest decay point of the 
decay straight L — d + s, i.e., the full flux (integrated over the angle) is seen by the 
detector from any decay point in the straight. 

Far detector limit In this case, the beam diameter given by the opening angle in Eq. (6) 
applied to Eq. (7) is of the order of the detector diameter at the nearest decay point 
L — d, i.e., the far distance approximation in Eq. (8) (first condition) is approximately 
fulfilled for any decay point in the straight. This implies that the beam spectrum will 
be similar to that of a far detector. 

Our near detector parameters are shown in Table 2. Note that the active volumes are much 
smaller than the actual detectors, because the secondary muons have to be stopped. There 
we define a hypothetical ND1 as a detector operating in the near detector limit, and a ND2 
as a detector operating close to the far detector limit. For example, we find from Eq. (7) 
using 9 in Table 1 (A = 0.9) that the beam diameter is about 17 m for L = s + d = 680 m, 
which explains the large diameter of ND1 to catch the whole flux. ND3 is an intermediate 
case between the near and far limits, as we will demonstrate later. The size of ND2 and 
ND3 is similar to conventional near detectors, such as SciBOONE, MINERz/A, NOMAD or 
the MINOS near detector. ND4 is a smaller version of ND2 with the same ratio between 
detector diameter and distance d. If the source was a point source (i.e., the straight would 
be a point), the event rate would be exactly the same as in ND2. Finally, ND5 is an 
OPERA-sized near detector, which we will only use for non-standard physics tests. For a 
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low energy, no storage ring is specified yet. Therefore, we assume the same geometry as in 
Fig. 3 (although a smaller storage ring could be sufficient). Since the beam is much wider 
because of the smaller boost factor (c/., Table 1), ND3 will already perform similar to a far 
detector. Therefore, we only use ND3 for this option. Of course, one can always up-scale 
or down-scale the discussed detectors (to re-scale the event rate), or change the distance 
to the source. However, from this qualitative discussion, we have covered all relevant cases 
allowed by the source geometry. In fact, ND4 may come closest to a realistic near detector. 



3.2 Near detector fluxes 



In order to compute the fluxes for the near detectors, let us first of look at the number of 
muon neutrino events produced in a near detector. From the very first principles, it is given 
by 

dN _ a dN Beam 1 f d 2 <$> M Det 

dE-A^^ NDct - a A^ t J dEdA dA ^' (9) 

where A^eam is the number of neutrinos in the beam (within the detector), Afoet is the 
number of target nucleons in the detector, a is the cross section per nucleon, and is the 
nucleon mass. The flux $ is, for a cylindrical detector, related to the decay rate T by 

d^ _n» d 2 v i (1Q) 



dEdA T dEd cos 9 2nL 2 

with the number of useful muon decays, T total decay rate, and L the baseline. Applying 
the far distance and point source approximations Eq. (8) in Eqs. (9) and (10), we have for 
the "point source" (PS) event rate 



%3 

/ 



° ~a — / "7771 77 1 6»=o sin6> d6 — ~ — - a 



dE r " A Dct J dEdcosO'" u ~ ™" m N r " dEdcos6' d m N 2ttL 2 



dN PS n M 1 f d 2 T . „,„M Det rv d 2 T M Dct 1 

~ — IT / I /i a sin H rff) ~ — — it 1 n „ 

v 2ttL 2 ' 

(11) 

where S mffi ~ D/{2L) is related to the detector diameter D. Therefore, only the on-axis flux 
is needed, the event rate increases linearly with detector mass, and drops as 1/L 2 . Such a 
flux is typically used for long-baseline experiment simulations, such as in GLoBES [52,53]. 

In our case, we cannot use these approximations for the near detectors. Therefore, we 
proceed in two steps: First, we fix the production point and integrate over the surface area 
of the detector. And second, we integrate over the decay straight. For a fixed production 
point in a distance L from the detector, we obtain similar to Eq. (11) 

D 
2L 

dN n„ 1 f d 2 T . „ ,„ Mtw dN™ A^ dN P 



^ 1 f d 2 Y . a , a M Bct dN PS A cR dN PS 

— 77^°"-; — ~n=r, -smOdO = — - — = — e(E,L) (12) 

T A Dct J dEd cos9 m N dE A Dct dE v 1 v ' 

o 

with the effective area A cf[ and the efficiency ratio e(E, L) determined by 



icff 



D 

2:TU ~ I a ^ 6de and e(E,L) = 4^, (13) 

J dEdcosB A Det 

o 



dEd cos 
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respectively. Therefore, the event rate in an arbitrary near detector can be related to the 
point source event rate in Eq. (11) with the same L. Then the efficiency ratio e describes 
what fraction of the beam is captured compared to the on-axis flux. Note that the event rate 
is highest on-axis for high neutrino energies and off-axis for low neutrino energies, which 
means that e can be larger than one for low neutrino energies, because the detector then 
captures some of the low energy part which is not present in the point source approximation 
(c/., Fig. 1: for the off-axis spectra, the low energy part is enhanced). From Eq. (13), one 
can read off that A cS approaches A Dct , or e goes to unity, if the far distance approximation 
in Eq. (8) is applied. 

As the next step, we average Eq. (12) over the decay straight assuming that the probability 
for muon decay is the same everywhere in the straight (cf., Fig. 3): 

d+s d+s 

«w- _ W ™ dL = i / ^! £M)i[ . (14) 



dE s J dE " s J dE 

d d 

From Eq. (11), we find that dNp ^' E "> qq \/ L 2 . Therefore, we can pull out dNvs ^> E ) f rom the 
integral in order to obtain 

dN wg _ dN PS (L cfi ,E)L 2 c{i d fe(L,E) _ dN PS (L cS ,E) 



dE dE s 

d 



j^ dL ^J^) m (15) 



with the average efficiency ratio 



s 



d+s 

L 2 

d 



HE) = & / E -^dL. (16) 



If we choose L eS = sj d(d + s) as the geometric mean between the farest and nearest point 
baselines of the production straight, we have that e(E) approaches unity if e(L, E) = 
1 or i > s (far distance or point source approximation in Eq. (8)). The meaning of 
Eq. (15) is the following: Implementing a near detector, we can use the same flux as for 
a conventional far detector, but with an effective baseline L e ff = a/ d(d + s), multiplied by 
the energy- dependent efficiency ratio e(E) to be computed from the geometry of the source 
and detector. Therefore, for instance, for an implementation in the GLoBES software, only 
the effective baseline and e(E) is needed. 

We show in Fig. 4 the efficiency ratio e(E) for the nearest point of the decay straight to 
the detector, the furthest point, and the average i(E) for ND1, ND3, and ND4 and for 
(the efficiencies are flavor dependent). Obviously, the efficiency ratios decrease with 
energy, since for large energies the beam is in most cases smaller than the detector. For 
low energies, however, a part of the off-axis flux can be captured, leading to an increase of 
events compared to the on-axis case. The discontinuities in Fig. 4 come from Eq. (3): If 
the beam, which has an energy dependent spread, becomes smaller than the detector, the 
efficiency ratio (or effective area) strongly decreases. From Fig. 4, we can read off that ND4 
(and, similarly, ND2, which we do not show) performs similar to the far detector limit, since 
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Figure 4: The efficiency ratios as a function of E v for ND1, ND3, and ND4. The curves for the farest 
point (efp) an d the nearest point (efp) to the near detector, as well as for the averaged efficiency e are 
shown. The horizontal dotted lines express the far detector limit (same spectrum as on- axis spectrum), and 
the dashed curves the near detector limit (whole flux captured). 





Neutrinos (z/^) 


Antineutrinos (z/ M ) 


ND1 


1.71 x 10 10 


8.82 x 10 9 


ND2 


3.69 x 10 9 


1.93 x 10 9 


ND3 


1.39 x 10 10 


7.21 x 10 9 


ND4 


9.36 x 10 8 


4.89 x 10 8 



Table 3: Total event rates for ND1 to 4 for a high neutrino factory with = 25 GeV in ten years (see 
definition in Sec. 4). 

the efficiency ratios are close to one (even in the high energy part, at least within the energy 
resolution of the detector). ND1 has efficiency ratios strongly decreasing with energy, which 
means that the low energy part of the spectrum becomes enhanced. It corresponds to the 
near detector limit (whole beam captured, c/., dashed curves), as anticipated. ND3 is an 
intermediate case: The farthest part of the decay straight resembles to the far detector 
limit, the nearest part the near detector limit. 

In Fig. 5 the double differential decay rates similar to Fig. 1 are shown, but multiplied 
with the corresponding efficiency factors e. Obviously, compared to the far detector case, 
the fluxes peak at lower energies, and the low energy part may be enhanced. Again, ND4 
performs close to a far detector, ND1 peaks at lower energies. The other two near detectors 
are intermediate cases. The event rates for ND1 to 4 are given in Table 3. Obviously, 
these rates are much higher than the ones in the far detector (a few hundred thousand in 
a distance L = 4000 km), which means that the detectors can in practice be built smaller 
(unless needed for other measurements, such as for particular cross section measurements). 
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Figure 5: The spectra (differential decay rates multiplied with e) for ND1, ND4, and the far detector 
(FD) limit, where the solid curves refer to the flux and the dashed curves to the v e flux. 

4 Systematics treatment and simulation methods 

Here we shortly describe the IDS-NF baseline setup including its systematics. Then we 
summarize our new systematics treatment, where details can be found in Appendix A. 
Finally, we describe our simulation methods. 

4.1 IDS-NF baseline setup 

The currently discussed setup for a high energy neutrino factory is the IDS-NF baseline 
setup [44]. It consists of two baselines L\ ~ 4 000 km and L2 — 7500 km, operated by 
two racetrack-shaped storage rings simultaneously with both polarities and \T stored, 
circulating in different directions). The geometry of a storage ring is shown in Fig. 3. There 
are no near detector specifications yet. We define the polarities of the muons stored in the 
rings as: 

+ : \i~ — > e~ + v e + ff, (17) 
- : /i + -> e + + v e + P„ (18) 

The following oscillation channels are used (with the corresponding muon polarities): 

Up appearance (— ) : v e — > (19) 

appearance (+) : v e — > (20) 

disappearance (+) : — > (21) 

disappearance (— ) : — > (22) 

For the backgrounds, neutral currents are included for all channels, and mis-identified 
charged current events (from the disappearance channels) are included for the appearance 
channels. The background levels are about 10~ 4 , depending on the energy (cf., Ref. [44] for 
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details). Since there are two racetrack-shaped storage rings SI and S2 targeted towards 
two far detectors, there are altogether eight oscillation channels. In the IDS-NF baseline 
setup 1.0, the systematics treatment is rather straightforward. For each channel and base- 
line, an overall normalization error is included, which is 2.5% for the signal rates, and 20% 
for the background rates. The normalization errors are uncorrelated between signal and 
background, among different channels, detectors, and polarities, but fully correlated among 
all bins. The energy resolution is assumed to be AE [GeV] = 0.55 \J~E~ [GeV]. 

4.2 Refined systematics treatment 

We refine the IDS-NF systematics treatment by introducing the following systematical er- 
rors, focusing on the cross section measurement as main purpose of the near detectors: 

Flux normalization errors, fully uncorrelated among the different polarities +, — and 
storage rings SI, S2, but fully correlated among all bins and all channels operated 
with the same beam. For example, the flux normalizations in Eq. (19) and Eq. (22) are 
treated fully correlated, since they come from the same fi + decay straight. Therefore, 
there are four independent errors (two storage rings times two polarities). We assume 
that the fluxes are known up to 2.5% without knowledge from the near detectors, 
which is a very conservative starting hypothesis from the current IDS-NF systematics 
assumptions. However, we also test the impact of an improved flux knowledge of 0.1%, 
which may be obtained by using various beam monitoring devices [42]. 

Cross section errors for the inclusive charged current cross sections, fully correlated among 
all signal and background channels measuring or v^, but fully uncorrelated among 
all bins. 3 We use bin widths of 1 GeV from 1 to 10 GeV, of 2 GeV from 10 to 20 GeV, 
and of 2.5 GeV from 20 to 25 GeV, i.e., 16 bins in total. This means that we at 
least crudely follow the energy resolution of AE [GeV] = 0.55^/E [GeV]. Since we use 
16 bins and and cross sections, there are 32 uncorrelated errors. We adopt a 
conservative point of view and assumed that the cross section for each bin is externally 
known to about 30% (see, e.g., summary plot in Ref. [54]). 

Background normalization errors, fully correlated among all bins, but fully uncorrelated 
among all channels, polarities, and detectors. 4 In total, there are two (polarities) 
times two (baselines) times two (channels), making eight background errors in the far 
detectors. We assume a 20% error each, just as in the IDS-NF baseline setup. 

Compared to the IDS-NF systematics, the different signal errors are correlated in a particular 
way. For instance, the cross section errors at the two far detectors are fully correlated, 

3 These errors are more accurately called "shape errors", because they are practically introduced after 
energy smearing (whereas cross section errors in principle enter before energy smearing). If the binning 
reflects the energy smearing, the cross section errors directly correspond to the shape errors. 

4 In principle, some knowledge on the backgrounds can be extrapolated from the near to the far detectors. 
However, the main backgrounds (charge identification and neutral current events) should rather depend on 
the charge identification and background discrimination properties of the far detector than the near-far 
extrapolation. Therefore, we use this conservative approach. 
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which will turn out to have interesting effects. For the efficiencies, however, we choose 
the IDS-NF numbers. Note that there might be some uncertainties coming from these 
efficiencies (such as a fiducial volume error), which we assume to be small and which we 
do not consider separately because they are not directly relevant for the near detector 
discussion. In addition, we neglect energy calibration errors. 

For our near detectors ND1 to 4, we introduce the same systematics for the backgrounds 
as for the far detectors, i.e., an uncorrelated 20% error for each polarity. For the sake 
of simplicity, our near detectors have the same energy resolution and binning as the far 
detector. However, as one important difference, we do not impose charge identification, 
since only the muon flavor is measured. Therefore, we assume 90% MINOS-like efficiency, 
starting at a threshold of 1 GeV. Only NC backgrounds are considered, at the level of the 
far detectors. For details of the systematics treatment, see Appendix A. 



4.3 Simulation techniques 



For the experiment simulation, we use the GLoBES software [52,53] with user-defined sys- 
tematics. For the oscillation parameters, we use (see, e.g., Refs. [1,55]) sin 2 # 12 = 0.3, 
sin 2 023 = 0.5, Am 2 x = 8.0 • 10~ 5 eV 2 , Am^ = 2.5 • 10~ 3 eV 2 , and a normal mass hierarchy, 
unless specified otherwise. We impose external errors on Am^ and 6*i 2 of 4% each, and we 
include a 2% matter density uncertainty [56,57]. 

The fluxes of the near detectors are included in GLoBES as user-defined fluxes. The user- 
defined fluxes are obtained from the product between the on-axis (normalized) rate Eq. (11) 
and the energy-dependent efficiencies in Eq. (16). Note that for the baseline, L e g has to be 
used. The event rate in GLoBES, similar to the one described in Eq. (11), is scaled auto- 
matically with the factors n M , a, and M^et- 5 In order to obtain the rest of the normalization 
in proper units, the normalization factor is determined as in App. C of the GLoBES manual 
in excellent agreement with the built-in fluxes: 

©norm = ^— — — x — — 5- x 1 [kton] x 10 9 [g/kton] x — l — x 10~ 38 [cm 2 ] 

27rx(l[km]) 2 10 10 [cm 2 /km 2 ] 1 J m N [g\ 1 J 

= 0.9581 x 10~ 16 . (23) 

Here is the average nucleon mass of the target material, i.e., l/(mjv[g]) is the number 
of target nuclei in one gram of target material. The first factor comes from the factor 
l/(27rL 2 ) in Eq. (11), i.e., the conversion from flux per angle to flux per surface area (the 
flux in GLoBES has to be specified in a distance of 1 km as flux per surface area). 



5 Standard oscillation physics with a high energy neutrino factory 



In this section, we perform a qualitative discussion of the impact of our systematics on 
different measurements, and we show several quantitative examples. 

5 For the sake of simplicity, we have not included time in the discussion, but only discuss the time- 
integrated event rates. 
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Figure 6: Simplified systematics for one storage ring. Here JQ refers to the cross section, to the 
cross section, and F\ and F2 to the flux normalizations. The index i marks the bin-dependence of the 
cross sections, whereas the flux normalizations are assumed to be bin-independent. 

5.1 Qualitative discussion 

For the following discussion, it is useful to have a simplified picture of systematics, which 
we illustrate in Fig. 6 for one storage ring, i.e., one baseline. Here Xi refers to the cross 
section, Xi to the cross section, and F\ and F 2 to the flux normalizations for n + and 
H~ stored, respectively. The index i illustrates the bin-dependence of the cross sections, 
whereas the flux normalizations are assumed to be bin-independent. The figure shows what 
products of cross sections and fluxes are measured in the different detectors. 

From the physics point of view, one may distinguish three different classes of measurements, 
which we have also tested quantitatively: 

Measurement of the leading atmospheric parameters. The main sensitivity for 
the leading atmospheric parameters Am§ x and 823 comes from the disappearance channel. 
From Fig. 6, we read off that only the product Fi Xi (or F 2 Xi for the different polarity) is 
measured in both the near and far detectors. Therefore, it is not necessary to disentangle 
Fx from Xi, and the product is precisely determined by the near detectors in each energy 
bin. For a large enough near detectors, this measurement is hence limited by the statistics 
of the far detector. However, without near detectors, the shape error coming from the cross 
section uncertainties will highly affect the shape- sensitive disappearance measurement. 

Measurements at the sensitivity limit of the appearance channel. At the sensitivity 
limit of the appearance channel, such as for the 6x3 and mass hierarchy measurements, the 
measurement will be limited by the backgrounds and their uncertainties. Therefore, even 
large uncertainties in the cross sections or fluxes do not harm the measurements, and the near 
detectors will not be required within our systematics treatment (which we have explicitely 
tested). However, it could help to better understand the backgrounds, a subject, which is 
beyond the scope of this study. 

CP violation measurement with the appearance channel. In this case, one has 

to distinguish the (moderately) small 8x3 and large #13 limits. In the small 8x3 case, the 
measurement is background dominated, where "background" refers to both systematics and 
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the CP-invariant contribution to the appearance probability (see, e.g., Eq. (31) in Ref. [58]). 
This contribution depends, for instance, on 6*23. Since the leading atmospheric parameters 
cannot be determined precisely without near detectors (see above), there will be some use of 
the near detectors no matter if the flux is known or not, i.e., the cross section uncertainties 
will be translated into a background uncertainty. 

For large enough statistics in the appearance channel, the knowledge on the signal normal- 
ization (such as on cross sections and fluxes) becomes an important factor. In this case, 
the products F\Xi and F2 Xi will be determined by the near detectors (c/., Fig. 6). How- 
ever, at the far detector, the different products F 1 Xi and F 2 Xi are needed to extract the 
information. Therefore, the ratio between the neutrino and antineutrino events, which is 
directly related to the CP violation effect, can be written as 

K (F\Xj) P{v e ^v,) (FjJQ X? P[v e 

°" {F 2 Xi) P(u e - Vp) (F 2 X,) X? P(P e 

FD measured ND measured 

where the first factor in the last product is measured by the near detectors. This means 
that only the product Xf/Xf x P(u e — > v^)/P(i) e — > P^) can be extracted, which leads to 
a different CP violation interpretation depending on the cross sections. In this case, there 
is an anti-correlation between fluxes and cross sections determined by F\ Xi — C± — const. 
and F2 Xi = C2 = const, in the near detectors independent of the energy bin %, and the cross 
sections can only be extracted from these relationships to the degree the fluxes are known. 
In the following, we assume that the fluxes can be monitored at the level of 0.1% if the 
near detectors are present, i.e., this anti-correlation can be broken. For instance, the 
event rates from purely leptonic elastic scattering interactions, for which the cross sections 
can be exactly calculated, should be sufficiently high to provide this precision [42]. 

5.2 Measurement of the atmospheric parameters 

As already discussed above, the near detectors are very important for the atmospheric 
parameter measurements, which are very sensitive to spectral effects. The impact of the 
near detectors is illustrated in Fig. 7 for a single baseline neutrino factory (L = 4000 km). 
The unfilled contours represent the far detector only, whereas the filled contours represent 
the far-near combination. In the right panel, the effect of the near detectors is very large 
for maximal mixing. Note that the unfilled contours depend on the knowledge of the cross 
sections, whereas the filled contours are limited by the statistics in the far detector. In the 
left panel, a non-maximal value of sin 2 #23 is chosen. In this case, the effect is less dramatic 
than in the right panel for maximal mixing, but still substantial. Especially, the octant 
degeneracy can be excluded with the near detectors at a high confidence level (if sin 2 2#i 3 
is large enough). 

The results hardly depend on which of ND1 to ND4 is chosen, because they all have sufficient 
statistics in all bins. This is illustrated in Fig. 8, where the Am^ allowed range is shown 
as a function of a luminosity scaling factor rescaling the near detector masses x operation 
time for two of the near detectors (ND3 and ND4). Comparing the large number of events 
in our near detectors in Table 3 with the few hundred thousand events in the far detector 
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Figure 7: The sin 2 823- Am^ allowed region for a high energy neutrino factory at the L = 4000 km 
baseline only; la, 2a, 3cr CL (2 d.o.f.), best-fit points marked by diamonds. The filled contours correspond 
to our near detector-far detector simulation, whereas the unfilled contours represent the far detector only. 
Normal hierarchy only, i.e., no sign-degenerate solution shown. In the left panel, sin 2 2^13 = 0.08 and 
<5cp = 0, in the right panel, sin 2 2#i3 = 0. 

(disappearance channel), it is not surprising that even much smaller near detectors would 
do. Basically, the luminosity starts to become important at the point when the near detector 
rates are of the order of the far detector rate. From Fig. 8 and Table 2, we can read off 
at the example of ND4 that, in principle, near detectors with a fiducial mass below 1 kg 
are sufficient for the standard oscillation parameter measurements. In this case, the near 
detectors can be easily operated in the far distance limit, i.e., the first part of Eq. (8) 
applies. Because there is no significant difference among the different near detectors, we 
only refer to the "near detectors" in the following. 

Very interestingly, if the neutrino factory is operated with two baselines at L\ = 4000 km and 
L2 = 7500 km, there is no significant effect of the near detectors - despite the fact, that two 
additional flux normalization errors for the second storage ring are introduced. Obviously, 
the spectral errors can be resolved by measuring the same cross sections in two different far 
detectors similar to a near-far combination, and energy- independent normalization errors 
are of secondary importance. 

5.3 CP violation measurement 

We show in Fig. 9 the CP violation discovery reach as a function of true sin 2 2# 13 and the 
fraction of (true) 6cp for one far detector (left) and two far detectors (right). The dashed 
curves represent the IDS-NF baseline setup with the corresponding systematics treatment, 
the solid lower curves our systematics treatment without near detectors and flux uncertain- 
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Figure 8: The allowed range for Amf^ as a function of a luminosity scaling factor rescaling the near 
detector masses x operation time (for ND3 and ND4). The contours are shown for the la, 2cr, 3cr CL (1 
d.o.f.). For L\ = 4000km and normal hierarchy only (i.e., no sign-degenerate solution shown), for maximal 
mixing and sin 2 2(?i3 = 0. 



ties (Jpiux = 2.5%, and the solid (thick) upper curves our systematics treatment with near 
detectors. In the left plot, we include a curve with a better known flux, whereas in the 
right plot, this curve coincides with the thick curve. As first important observation, our 
systematics treatment (including near detectors) leads to a better sensitivity that the IDS- 
NF systematics, in spite of our conservative choices for the systematical uncertainties. The 
main difference between these two systematics treatments is that we assume the systemat- 
ical errors to be correlated among different detectors, whereas the IDS-NF setup assumes 
fully uncorrelated errors. The practical realization will likely be somewhere in the middle 
(between the dashed and thick curves), since, for instance, cross sections are fully correlated 
among all detectors measuring the same flavor and polarity, whereas there may be other 
normalization and spectral errors which depend on the detector (such as fiducial volume 
errors leading to errors in the detector-dependent efficiencies). Similar to the atmospheric 
parameter measurement, the near detectors are important for one baseline only, whereas a 
two-baseline neutrino factory hardly benefits from the near detectors. 

For one baseline (left plot), an improved flux knowledge only helps for large # 13 in the 
absence of the near detectors, because the correlation described in Sec. 5.1 can be resolved 
and only the uncertainty of the atmospheric parameter measurements remains. We have 
also tested a better matter density knowledge of 0.5% (compared to 2%), which slightly 
improves the measurement in all cases for large #13. In the case of two baselines with near 
detectors in the shorter baseline storage ring (thick curve in right plot), we have found no 
improvement from better known fluxes. Since the cross sections can in this case already 
be extracted by the two near detectors in the first storage ring, this implies that putting 
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Figure 9: CP violation discovery reach as a function of true sin 2 2#i3 and the fraction of (true) Sep for 
one far detector (left) and two far detectors (right); 3<r CL. 

additional near detectors in front of the second storage ring, such as for flux monitoring, 
may not be required from the physics point of view. From the figure, neither a precise flux 
monitoring nor near detectors are mandatory for a successful CP violation measurement in 
a two baseline neutrino factory, if the far detectors are well enough understood. 

6 Standard oscillation physics with a low energy neutrino factory 

Here we briefly discuss the use of near detectors for a low energy neutrino factory. For the 
low energy neutrino factory, there is not yet any baseline specification within the IDS-NF. 
Therefore, for our reference setup, we follow Ref. [51], using = 4.12 GeV and their "high 
statistics" scenario with 10 23 kt yr exposure, corresponding to 10 years of data taking with 
5-10 20 useful muon decays per polarity and year, and a fiducial detector mass times efficiency 
of 20 kt. We use the baseline L = 1290 km, which corresponds to Fermilab-Homestake. 
The detector proposed in this reference is a magnetized TASD (Totally Active Scintillating 
Detector) with a threshold of about 500 MeV, an energy resolution AE [GeV] = 0.3 E [GeV] 
(conservative estimate), and an estimated detection efficiency of 73%. For the binning, we 
use nine equi-distant bins between 500 MeV and 5 GeV. The channels are the same as 
used for the high energy version, but the background level is conservatively estimated to be 
10 -3 [51]. Note that we include two types of backgrounds for the appearance channels, one 
which scales with the disappearance rates (such as from charge mis-identification), and one 
which scales with the un-oscillated spectrum (such as from neutral current events), both at 
the level of 10~ 3 . As in Ref. [51], we take the systematical errors to be 2% for all signal and 
background errors in the reference setup. We can reproduce the event rates in Ref. [51] very 



18 




Figure 10: The sin 2 ^23-^31 allowed region for a low energy neutrino factory; la, 2a, 3a CL (2 
d.o.f.), best-fit points marked by diamonds. The filled contours correspond to our near detector-far detector 
simulation, whereas the unfilled contours represent the far detector only. Normal hierarchy only, i.e., 
no sign-degenerate solution shown. In the left panel, sin 2 2613 = 0.08 and Sqp — 0, in the right panel, 
sin 2 20i3 = 0. 

well, as well as the fits in the ^-Am^-planes and #i3-<5cp-planes for their set of parameters. 
Compared to Ref. [51], we do not find any sensitivity to the mass hierarchy for 6*13 = 0. 

For the near detector implementation, we follow the same procedure as for the high energy 
neutrino factory using ND3. However, we only have one storage ring and one far detector, 
and nine (equi-distant) bins, making a total of 18 cross section errors (instead of 32). The 
external knowledge of the fluxes ofi ux is assumed to be 2% (if not stated otherwise), and 
the backgrounds are assumed to be known at the same level [51]. 

In Fig. 10, we show the results for the leading atmospheric parameters for a deviation from 
maximal mixing and large #13 (left), and maximal mixing (right). The unfilled contours 
correspond to no near detector, and the filled contours include the near detectors. As 
for the high energy neutrino factory with one baseline only, the near detectors are very 
important for the atmospheric parameter measurement. Of course, the results without near 
detectors will depend on the assumptions on the cross sections. However, one can easily 
see the benefits of the near detectors, especially since there is typically no second baseline 
planned for the low energy neutrino factory. 

In Fig. 11, we show the main results for CPV (where the effect is most dramatic). Obviously, 
the setup with near detectors corresponds very much to the reference setup [51], except for 
large 9 13 , where the reference setup is slightly better. Without near detectors, the result 
is significantly worse, which means that the near detectors are mandatory for a low energy 
neutrino factory of all #13 to better understand the cross sections. In either case (with 
or without near detector), a better known flux (0.1%) helps somewhat for large sin 2 2#i3, 
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Figure 1 1 : CP violation discovery reach as a function of true sin 2 29^ and the fraction of (true) <5cp for 
a low energy neutrino factory; 3tr CL. The reference setup is taken from Ref. [51]. 

whereas the matter density uncertainty is of secondary importance for this setup. As for the 
high energy neutrino factory with only one far detector, a better known flux yields almost 
no improvement for small #13. 

7 Non-standard interactions with a high energy neutrino factory 

Non-standard interactions (NSI) [59-63] are effects of physics beyond the Standard Model 
(SM) which can be described by four fermion vertices below the electroweak symmetry 
breaking scale relative to the SM interactions: 

6J? cS = 2V2G F (e L /% 7 (i/VP^) (£ 6 YP l/ r£ 7 ) • (25) 

Here Pl and Pr are the the left- and right-handed (chiral) projection operators, respectively. 
The operator in Eq. (25) may come from effective non-renormalizable d = 6 [64-66], d = 

8 [67,68], etc., operators using the SM field content (including the Higgs). The e from 
the dimension d operator is suppressed by e oc (w/v / 2A) d_4 with respect to the SM Higgs 
vacuum expectation value v by the new physics scale A, which comes from integrating out 
the heavier (new) fields. 

These NSI can affect the neutrino propagation by coherent forward scattering in Earth 
matter if the two charged leptons in Eq. (25) are electrons, which is only sensitive to the 
vector component as 

= with eZ L/R = {e L/R )%. (26) 
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In addition to the propagation in matter, the production or detection processes can be 
affected by NSI. For the specific case of a neutrino factory and considering just purely 
leptonic NSI, dominantly effects at the source are relevant, since the (initial) detection 
interactions involve quarks. They are customarily parameterized in terms of e s a/3 , which 
describes an effective source state (for small e s ) as [17,62,69] 

K) = K) + E <#\ V f>) With <U = ^)% OT <3=(eX" ( 27 ) 

f3=e,/j,,T 

In this case, the muon decay rate could be modified by the NSI interaction in Eq. (25), with 
the largest effect resulting from the coherent contribution to the state at the source [17,22]. 
It appears as an admixture of a given flavor v a with all other flavors, encoded by v$ in 
Eq. (27). The dominating effect comes from the left-handed component, since the right- 
handed component will be helicity suppressed at production. 

Near detectors may be very relevant for the extraction of the source NSI, since these lead 
to an unambiguous "zero-distance" (L <C L osc ) effect proportional to |e s | 2 . However, within 
our systematics framework, how should one disentangle source NSI modifying the electron 
or muon flavor contributions (e s ee , e s , e* e , and e* ) from the unknown fluxes and cross 
sections? For example, a larger than anticipated muon neutrino rate at zero distance may 
come from underestimated cross sections, the neutrino flux normalization, or e s , i.e., it 
will be very difficult to disentangle these contributions, especially if more of the e's are 
allowed simultaneously. For e s efl , however, near detectors with charge identification could 
help (because there will be of the opposite polarity in the beam), and for e* e near detectors 
with electron flavor measurement and charge identification. Since our near detectors are 
more primitive than that, we do not consider these measurements. However, note that 
for some new physics measurements, such as electron neutrino disappearance into sterile 
neutrinos [70], these additional channels are important. 

Because there are no tau neutrinos in the beam, the most interesting option may be to use 
near detectors to measure v T appearance, see, e.g., Refs. [37,71]. Therefore, we use the 
high energy neutrino factory with ND3 and two (symmetrically operated) additional near 
detectors at L — 1 km for v T detection. In the absence of detailed studies for alternatives at 
a neutrino factory, we adopt OPERA-like properties with 2 kt detector masses, as simulated 
in Ref. [10] and described in Ref. [8] for the v e — > v T channel at the neutrino factory. These 
detectors correspond to ND5 (see also Fig. 3). Since there is no study for — > v T yet, we 
follow the pragmatic approach in Ref. [72] (see also Ref. [73] for a similar approach) and 
use five times the signal and five times the backgrounds from Ref. [8], because no charge 
identification is required, and therefore all tau decay modes (the hadronic ones as well) may 
be exploited at the price of higher backgrounds. We assume that similar efficiencies and 
backgrounds can be used for both the neutrino and antineutrino channels, and we add all 
events from u e — > v T and — > v T (same for the other polarity). For the detector geometry, 
we use OPERA-like measures, meaning that the detectors will perform as ND5 from the 
flux point of view. We show the differential decay rates corrected by the efficiency factors 
for these detectors in Fig. 12. Compared to the far distance limit, which is also shown 
for comparison, there is already a substantial modification of the fluxes. Since we do not 
measure the v T cross sections explicitely, we assume increased 30% signal normalization 
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Figure 12: Spectra (differential decay rates multiplied by the efficiency factor £) for the and v e fluxes 
of the OPERA-sized ND5 (d = 1 km). For comparison, the spectra of the far distance (FD) limit are shown. 

errors. However, at the sensitivity limit for e s , the signal normalization errors will be of 
secondary importance, meaning that the cross sections do not need to be exactly known 
for this application. Note that this toy implementation of a u T detector may be in fact too 
primitive in practice, it should only serve as an example. For instance, the scanning load 
in an OPERA-like emulsion cloud chamber might be extremely high, which means that 
a different detector technology with at least similar efficiencies and vertex reconstruction 
abilities might be desirable {e.g., liquid argon). Maybe the and v T near detection can 
even be done at the same place with the same detectors. 

We give the sensitivities for the source NSI |eg T | and |e* |, and the matter NSI |e™|, and 
|e^.| in Table 4 with and without v T near detectors. As expected, the source NSI improve 
significantly in the presence of the near v T detectors - especially |e* |. Compared, for 
instance to the NOMAD results [74] of 0.086 for \e% T \ and 0.013 for |e* T | (90% CL), the 
results can be improved by about two orders of magnitude. 6 The bounds of the matter NSI 
|e™| and |e^.| compare very well to the ones in Ref. [76] even for our different systematics 
treatment. Note that the bound on |e^_| is rather weak if e™ is allowed to be complex. If, 
however, the NSI are assumed to originate from effective dimension six operators involving 
four lepton doublets, and charged lepton flavor violation is sufficiently suppressed, there is 
a general connection between source and matter NSI as [77] 

£r = -(#)'. (28) 

This means that the bounds on |e* | applies, under this condition, to \e™\ as well, and vice 
versa. In particular, the bound for \e™\ becomes quite strong in this case, even better than 

6 There are also other theoretical limits, such as from rare lepton decays [75]. However, the neutrino 
factory limits will be even stronger. 
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Table 4: Sensitivity to the NSI parameters without (left column) and with (right column) 
v T near detectors at the 90% confidence level. For the non-diagonal NSI, the phases have 
been marginalized over. 

the bound from lepton universality 1.9 10 -3 (90% CL) [78] in the presence of the v T near 
detectors (see last row in Table 4). 

An interesting application of the correlation in Eq. (28) is the test of non-standard CP 
violation (NSI-CPV). Such NSI-CPV has been studied in Refs. [17,72,79-81]. Since the 
zero-distance effect is proportional to |e s | 2 , there is no sensitivity to CP violation in our 
near detectors. However, if the source and matter NSI are correlated through Eq. (28), the 
near detectors can be used to constrain e* T , and the far detectors can be used to measure 
the NSI-CPV in e™ very well (as discussed in Ref. [80]). In summary, our configuration can 
identify such NSI-CPV for <j>™ = tt/2 (maximal CPV) in that case down to |e™ | ~ 0.0005 
(90% CL), whereas without v T near detectors it can only be measured down to 0.003. This 
is comparable to the result in Ref. [72], where the additional baseline L = 130 km was 
proposed. Without the correlation in Eq. (28), the CPV in e* T can be measured down to 
0.001 in the far detectors (for = — 7r/2), meaning that the additional contributions from 
e™ T are be counter-productive in the presence of the correlation (see Fig. lb in Ref. [22]), 
and the CPV in e™ down to 0.02 only, in consistency with Ref. [80]. 

8 Summary and conclusions 

We have discussed the impact of near detectors on standard oscillation physics at a high and 
low energy neutrino factory, and the possible use of near detectors for u T appearance. We 
have included the geometry of the source and the detectors in the event rate calculations, 
and we have defined qualitatively different near detectors. For this definition, it has turned 
out to be useful to consider the near detector limit (detector captures whole beam flux) and 
the far detector limit (detector observes similar spectrum to that of a far detector), as well 
as intermediate cases between these limits. We have demonstrated that a neutrino factory 
near detector can effectively be parameterized as a far detector with an energy-dependent 
effective area. In this case, the decay straight can be treated as a point source with an 
effective baseline which is the geometric mean between the shortest and longest distance to 
the detector. 

In order to test the impact of the near detectors on physics, we have used the IDS-NF 
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(International Design Study for the Neutrino Factory) baseline setup as a starting point. 
In the current baseline setup, the near detectors are not yet defined explicitely. Therefore, 
we have refined the systematics treatment including cross section errors fully uncorrelated 
among the energy bins, but fully correlated among all channels measuring muon neutrinos 
or antineutrinos. In addition, we have included flux normalization errors fully correlated 
among all energy bins, but uncorrelated among the different storage rings and muon polari- 
ties. As a conservative assumption, we have kept relatively large (uncorrelated) background 
uncertainties in all channels. We have then tested the impact of this systematics treatment 
and the different near detector types on standard oscillation physics, and we have considered 
u T appearance in near detectors for non-standard interactions. 

For the characteristics of the near detectors for standard oscillation measurements, we 
have found that two near detectors should be operated on-axis (with respect to the decay 
straights) if the muons and anti-muons circulate in different directions in a racetrack-shaped 
ring. These two detectors measure the products between the (uncorrelated) flux normal- 
izations and muon or muon anti-neutrino cross sections. For a two baseline operation of 
a high energy neutrino factory with two storage rings, we have demonstrated that near 
detectors in one storage ring might be sufficient. The near detectors should be at least as 
good as the far detectors in terms of efficiencies, energy resolution, etc., but the technical 
details should be of secondary importance. The near detectors can be very small for the 
systematics considered as long as the event rates are similar to that of the far detectors - a 
few kilograms are basically sufficient. In this case, the near detectors can be operated quite 
closely to the source in the "far distance limit", which means that the observed spectrum 
is similar to that of a far detector. For larger near detectors, the event rates are very high 
in all energy ranges so that a possibly different spectrum does not significantly affect the 
near-far extrapolation even in the most extreme cases. 

As far as the standard oscillation parameter measurements are concerned, the results de- 
pend on the number of far detectors used. If the neutrino factory has only one baseline, 
such as a high energy (E^ = 25 GeV) neutrino factory with L = 4000 km or a low energy 
(Efj, = 4.12 GeV) neutrino factory with L = 1290 km, the near detectors are mandatory for 
the leading atmospheric parameter measurements. For the #13 and mass hierarchy discovery 
reaches, we have not found any significant impact, because these measurements are back- 
ground limited. For the CP violation discovery reach, the near detectors are important, be- 
cause the unknown atmospheric parameters lead to intrinsic unknown backgrounds from the 
CP-even terms, which can be controlled by the near detectors. For small sin 2 2#i3 < 10~ 2 , 
this effect hardly depends on the flux knowledge, which means that flux monitoring is of 
secondary interest. For large sin 2 2# 13 , however, fluxes monitored to the level of 0.1% help 
to understand the measurement. In particular, it is ensured that the CP violating effect 
does not come from a wrongly assumed neutrino-antineutrino cross section ratio. 

For a two baseline high energy neutrino factory, such as the IDS-NF baseline setup with 
L\ = 4000 km combined with L2 = 7500 km, we have demonstrated that the considered 
systematical errors cancel even without near detectors and good flux monitoring. In such 
a setup, the same product of fluxes and cross sections is measured in both far detectors. 
Therefore, a two baseline neutrino factory should be very robust with respect to systematics, 
no matter how many near detectors are used. Of course, in staging scenarios where only one 
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Figure 13: #13, CP violation, and mass hierarchy discovery reaches (3tr) as a function of the fraction 
of Sep and (true) sin 2 2f5i 3 . The solid curves represent our near detector- far detector simulation for the 
high (left curves) and low energy (right curves) neutrino factories, whereas the dashed curves represent the 
IDS-NF [44] and Bross et al [51] reference setups. 

baseline is operated first, near detectors are still required, as well as they may be needed 
for different purposes. 

Because of the particular correlations of the cross section and flux normalization errors, 
our high energy neutrino factory performs slightly better than the current IDS-NF base- 
line version, where all signal and background normalization errors are fully uncorrelated. 
Therefore, systematics may be actually quite conservative for the IDS-NF baseline. Note 
that we have not even considered a background extrapolation from the near to far detector 
(depending on the properties of the near detectors), but instead adopted the conservative 
IDS-NF uncorrelated errors. However, we have not considered energy calibration errors, 
where the near detectors may also help, and detector-dependent systematical uncertainties, 
such as errors in the detection efficiencies, which may partially cancel if similar far detectors 
are used. Therefore, a more realistic systematics treatment may produce results somewhere 
in the middle. For the low energy neutrino factory, our systematics treatment reproduces 
the reference setup with effective near detector simulations very well. We summarize our 
main results in Fig. 13, where the # 13 , CP violation, and mass hierarchy discovery reaches 
are shown for both the IDS-NF and low energy neutrino factory. The solid curves represent 
our near detector-far detector simulation for the high (left curves) and low energy (right 
curves) neutrino factories, whereas the dashed curves represent the reference setups. This 
figure illustrates the impact of the different systematics, and it compares, for the first time, 
the low and high energy neutrino factory options for the same input assumptions. 

Finally, we have considered near detectors with OPERA-like properties for v T appearance 
to measure non-standard interactions (NSI). We have demonstrated the NOMAD limits for 
the zero-distance effect could be improved by about two orders of magnitude. In addition, 
we have discussed the theoretical possibility that source and matter NSI are connected at a 
generic level assuming dimension six effective operators. In this case, even the matter NSI 
parameter |e™ | becomes quite strongly limited, exceeding the bound from lepton universal- 
ity. Furthermore, CP violation from matter NSI may become measurable down to 0.0005 
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We conclude that two near detectors are mandatory for a successful neutrino factory oper- 
ation with one baseline, and a good enough flux monitoring will be useful for some physics 
measurements. However, a two baseline high energy neutrino factory may prove to be very 
robust with respect to systematical errors, even if the systematics goals in the initial one 
baseline operation phase cannot be achieved. In either case, near detectors will be manda- 
tory to constrain certain new physics effects. Since the statistics is typically very high, the 
detailed requirements for the near detectors will be most likely driven by flux monitoring 
and new physics measurements, rather than the cross section uncertainties. 
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A Details of the systematics treatment 

For the systematics treatment, we use the pull method. As an example, we describe here 
the systematics treatment of the high energy neutrino factory. The total x 2 is given by 

4 4 

X 2 = Yl Xfi,c + Yl x ^,c + X 2 ni + X%2 + Xpull , (29) 
C=l c=\ 

summed over the xjj c °f detector D and channel C. For the detectors, we have the far 
detector Fl at the 4 000 km baseline, the far detector F2 at the 7500 km baseline, the 
near detector Nl behind the fi + decay section of the storage ring SI to Fl, and the near 
detector N2 behind the yT decay section of the storage ring SI to Fl (c/., Fig. 3). For the 
far detector channels C, we have appearance (1), appearance (2), disappearance 
(3), and disappearance (4), whereas the near detectors only measure un-oscillated and 
Pfj,, respectively. 

Let the theoretical (fit) rate T % be the rate in the ith bin, which is composed of the signal 
rate S l and the background rate B l . Then we have for the different fit rates in the bins 
using the a's as the auxiliary parameters of the pull method to be marginalized over 7 

Tk,i = (1 + a&f + «S SU,! + (1 + 4g) *ki » ( 3 °) 
T k,2 = (1 + + 4"si) ^1,2 + (! + 4g) BU,2 , (31) 

T F1 ,s = (1 + + S F1 ,3 + (1 + 4g) ^i,3 > (32) 

7 Here we have simplified the background treatment a bit. In fact, we include the bin-dependent cross 
section uncertainties in the backgrounds from charge mis-identification. However, there is hardly any effect 
of this more refined treatment. 
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+ ^)S i F1A + (^ + a 4 BG )B i F1A , (33) 
+ Sk,i + (1 + 4g) B F2 ,i , (34) 

+ <h£) ^ 2 ,2 + (1 + 4g) ^2,2 , (35) 
(1 + <luf + ^2,3 + (1 + 4g) ^2,3 , (36) 

+ 4^c) Sk,4+ (1 + 4g)5K 4 , (37) 

+ 4^c) + (1 + «bg) , (38) 

+ a^: c )S i N2 + (l + a 1 B ° G )B i N2 . (39) 

The auxiliary parameters describe the following types of systematics errors: Op lu ' x 5fc are 
flux normalization errors for storage ring Sk and fj^ stored (fully correlated among all 
bins and detectors in the same beam), a^^ c represent the neutrino cross section errors 
(fully correlated among all channels measuring z/ M , but uncorrelated among the bins), axsec 
represent the antineutrino cross section errors (fully correlated among all channels measuring 
P M , but uncorrelated among the bins), and a BG the background errors (fully uncorrelated 
among all channels, but fully correlated among all bins). Altogether, there are 32 cross 
section errors, ten background normalization errors, and four flux normalization errors, i.e., 
46 systematical errors in total. 

The x 2 f° r each detector or channel is then obtained as 



rpi 

1 F1A ~ 


(1 + <f 


rpi 

1 F2,1 — 


(l + <f 


rpi 

_F2,2 — 


(1 + <iux 52 


rpi 

1 F2,3 — 
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rpi 

1 F2,4 — 


(1 + ^f 


rpi 

A N1 — 


(l + a&f 


rpi 

1 N2 — 
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V 1 D,C 



(40) 



i=i 

where the Poissonian x 2 , in which O is the observed/simulated rate and T is the theoreti- 
cal/fit rate, is used (near detectors analogously). In addition, the pull \ 2 is given by 



XPull 



* / fi+,Sk\ 2 2 / M -,5fc\ 2 16 / ^,i\ 2 16 / P M ,i\ 2 10 / k \ 1 

f a Fiux j _|_ y^ [ a Fiux j _|_ y^ [ a xscc j + y- a x SC c j _j_ y^ / %g \ 



i \ 0"Xscc / ~i \ CXsec / \0"BG/ 

1 \ / i=l \ / fe=l 

(41) 

For the high energy neutrino factory, we (conservatively) assume <tfi U x = 2.5%, oxsec = 30%, 
and o"bg = 20%, unless explicitely stated otherwise. 
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